Condensate deformation and quantum depletion of 
Bose-Einstein condensates in external potentials 



C A Muller 1 and C Gaul 2 

1 Centre for Quantum Technologies, National University of Singapore, Singapore 1 17543, 
Singapore 

2 GISC, Departamento de Fisica de Materiales, Universidad Complutense, E-28040 Madrid, 
Spain 

PACS numbers: 03.75.Hh Static properties of condensates; thermodynamical, statistical, and 
structural properties. 

Abstract. The one-body density matrix of weakly interacting, condensed bosons in 
external potentials is calculated using inhomogeneous Bogoliubov theory. We determine the 
condensate deformation caused by weak external potentials on the mean-field level. The 
momentum distribution of quantum fluctuations around the deformed ground state is obtained 
analytically, and finally the resulting quantum depletion is calculated. The depletion due to the 
external potential, or potential depletion for short, is a small correction to the homogeneous 
depletion, validating our inhomogeneous Bogoliubov theory. Analytical results are derived for 
weak lattices and spatially correlated random potentials, with simple, universal results in the 
Thomas-Fermi limit of very smooth potentials. 



1. Penrose-Onsager criterion and Bogoliubov theory for inhomogeneous condensates 

The phenomenon of Bose-Einstein condensation (BEC) underlies bosonic superfluidity and 
thus plays a pivotal role in condensed-matter physics. In recent years, the unique experimental 
possibilities offered by dilute ultracold atomic gases have triggered a renewed interest in 
Bose-Einstein condensates. These are loaded into external potentials of a large variety, 
ranging from simple harmonic traps to increasingly complicated optical lattices, all the way 
to random potentials [1—3]. The more complicated the confining potentials are, the greater is 
the challenge to tell the condensate from the excitations, both quantum and thermal. Yet, to be 
able to distinguish precisely between condensate and excitations is not only interesting from 
a conceptual point of view, it is also important in order to understand the physical properties 
of Bose-condensed gases. 

A criterion for BEC that applies to interacting as well as inhomogeneous systems was 
proposed in 1956 by Penrose and Onsager [ ] and has remained in vigor until today: BEC 
occurs whenever the one-body density matrix (OBDM) 

p(r,r') = (¥(r)4>(r')) (1) 

has (at least) one macroscopically occupied eigenmode. As stated very clearly by Penrose 
and Onsager in their original paper, only if the system is completely homogeneous (i.e. 
translation invariant under periodic boundary conditions), then condensation occurs into a 
single momentum component, namely the state with wave vector k — in the condensate 
rest frame. There is, certainly, numerical evidence that the k = component does have an 
important bearing on the superfluid property [5,6]. But unless the additional assumption of 



Bose-Einstein condensates in external potentials 



2 



spatial homogeneity is met, the zero-momentum occupation must not be used to determine 
the condensate fraction. Recently, Astrakharchik and Krutiksky [ ] have devised a quantum 
Monte Carlo scheme of computing the OBDM and condensate mode in external potentials. In 
the present article, we employ analytical Bogoliubov theory, applicable to weakly interacting 
condensates at low temperature, and calculate the condensate fraction and quantum depletion 
of inhomogeneous Bose gases. 

Bogoliubov theory describes quantum fluctuations around a mean-field condensate by 
splitting the quantum field 

*(r) = $(r) + 6$(r) (2) 

into a (large) mean-field condensate $(r) and (small) quantum fluctuations 5^(r). In the 
symmetry-breaking picture, the condensate $(r) = (^(r)) is the expectation value of the 
quantum field. By consequence, (5^(r)) = 0, and the OBDM (1) splits into the sum of a 
condensed and non-condensed contribution, 

p(r,r') = $*(r)$(r') + (5& (r)5$(r')) . (3) 

This form of the OBDM complies with the third version of the Penrose-Onsager criterion [ ], 
often quoted in a shortened manner. In full, this criterion reads as follows: If 

\p(r, r') - $*(r)$(r')| < wy(\r - r'\) (4) 

with a function •y(s) that is independent of the density n = N / L d and goes to zero at infinity, 
and if <l>(r) contains order N particles, then BEC occurs, and $(r) is a good approximation to 
the condensate wave function. An OBDM with the asymptotic form (4) is said to possess off- 
diagonal long-range order. Bogoliubov's ansatz holds whenever the condensed component is 
large, and then the non-condensed component of the OBDM (3) can be bounded as required 
by (4) [8]. 

For further analysis, it is convenient to consider the bulk-averaged OBDM 

p(s) = L~ d J d d r (#t( r )#( r + s )) = L- d c ik s n k (5) 

k 

that depends only on a single vector s. It is the (inverse) Fourier transform of the single- 
particle momentum distribution n k = (^j^fe)' where ^ k = L~ d l 2 J d d re~ lk ' r ^>(r) is the 
particle annihilator in momentum representation. Under the Bogoliubov ansatz (2), also the 
momentum distribution splits into a condensed and non-condensed component: 

n k = l^fel 2 + (£*fc£*k) = n ck + 5n k . (6) 

Starting from these well-known concepts, we investigate in the following the condensate 
deformation caused by weak external potentials, calculate the corresponding fluctuation 
momentum distribution, and finally determine the resulting quantum depletion. We apply our 
analytical theory to lattice potentials and spatially correlated random potentials. The paper 
is structured as follows. In Sec. 2, we recall the momentum distribution and OBDM of the 
mean-field condensate in presence of a weak external potential. In Sec. 3, we draw on the 
inhomogeneous Bogoliubov theory developed in [ ] and derive the momentum distribution 
of fluctuations in external potentials of arbitrary form. Notably, we discover a universal 
momentum distribution in the Thomas-Fermi regime of smoothly varying potentials. In 
Sec. 4, we compute the resulting quantum depletion. The depletion caused by the external 
potential is found to be a small correction proportional to the depletion of the homogeneous 
Bose condensate, thus validating the Bogoliubov theory of inhomogeneous condensates. 
Section 5 concludes. 
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Figure 1. (a) Condensate density n c (r)/n c , perturbatively calculated via (8) and (9), in a 
square lattice potential V(r) = V x cos(K x x) + V y cos(K y y) with V x = V y = 0.75gn c 
and K X S; = K y £ = 1. The condensate adapts a deformed configuation in the external 
potential by avoiding peaks and accumulating in wells, (b) Corresponding condensate 
momentum distribution n c k/N c , Eq. (11), showing 1/16 ~ 6% of the total population in 
the fc-components imprinted by the lattice. 



2. Condensate deformation 

2.1. Mean-field theory 

On the mean-field level, it is relatively simple to calculate the deformation of the condensate 
amplitude $(r) caused by an external potential V(r). One has to solve the stationary Gross- 
Pitaevskii (GP) equation [8] 

[(-fi 2 V 2 /2m) + V(r) + g\$(r)\ 2 ] $(r) = fJ,^{r) (7) 
at given chemical potential [i. Numerically, the solution $(r) can be computed rather 
efficiently by imaginary-time propagation of the time-dependent GP equation [10]. 



2.2. Condensate deformation by a weak potential 

When the external potential is weak, it is straightforward to solve the GP equation 
perturbatively to the desired order in powers of V [1 1, 12]: $ = + + $ (2) + . . .. 
In the following we work at fixed average condensate density n c = L~ d J d d r|$(r)| 2 and 
adjust the chemical potential accordingly [9]. In momentum representation, the homogeneous 
condensate $^ = Nc Skfi receives the lowest-order deformations 

$ L X) = -N^ /2 V k , (8) 
^ = 2an C +e° £ ^- q % [(1 - S k0 )e° q - gn c ] . (9) 

q 

The set of small parameters of this expansion are the reduced matrix elements 

1 — Si. o 

Vk = ■ 9 V k, (10) 

el + 2gn c 

with Vfc = L~ d J d d re~ lk ' r V(r), gn c — h 2 /(2m£ i 2 ) the mean-field interaction energy, and 
= h 2 k 2 /2m the bare kinetic energy. Consequently, the condensate momentum distribution 
up to order V 2 becomes 



n ck = |$ fe | 2 - N, 



(I - V 2 )6 k0 + \V k \ 



(11) 
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Most particles in this distribution still have zero momentum, but a small fraction 

V 2 = ^|Vfe| 2 «l, (12) 

k 

of them have been promoted to finite momenta by the weak external potential. 

As an illustration, Fig. 1 shows the condensate density n c (r) and its momentum 
distribution n c k, Eq. (1 1), in presence of a simple square lattice potential whose parameters 
are such that V<z — 1/16 = 6.25%. The external potential deforms the condensate, which 
tends to avoid potential peaks and accumulates in potential wells. The deformation is periodic 
in real space, Fig. 1(a). In fc-space, the lattice shifts some population to the lattice momenta, 
Fig. 1(b). To higher order in the external potential, also higher-order components would 
become visible in the momentum distribution (cf. Fig. 2 in [13]). 



2.3. Condensate deformation does not reduce mean-field condensate fraction 

Condensate deformation is clearly a mean-field effect. The condensed part of the OBDM 
follows by inserting (11) into (5): 

Pc (s) = n c (l-V 2 ) + n c J2^ k ' s \Vk\ 2 - (13) 

k 

By construction, p c (0) = n c is the total density of condensed particles, which is kept fixed 
by adjusting the chemical potential. Consequently, constant potential offsets have no effect, 
and thus Vq ~ as implied by (10). Since Vq = 0, the fluctuating part in (13) seems to tend 
to zero in the limit s — > oo, which would imply that it does not contribute to the off -diagonal 
long-range order. Then, the condensate fraction would only be (1 — V2), equal to the fraction 
of particles in the k = momentum eigenstate. However, this reasoning is erroneous, as 
becomes quite evident in the case of a general lattice potential 

d 

V(r) = Vj cosiKj ■ r). (14) 
3=1 

By consequence of the lattice periodicity, also the OBDM deformation is periodic in s, 

k j=l 

with finite amplitudes Vj = Vj/(e a K . + 2gn c ). But since a periodic deformation cannot be 
bounded by a function j(s) that tends to zero, it is the full momentum distribution (9) that 
describes the condensate, in agreement with the Penrose-Onsager criterion. And indeed, the 
(wrong) conclusion that only part of the field $(r) constitutes the condensate would contradict 
the mean-field ansatz used to calculate it in the first place. 

If one performs additional configuration averages, then one accesses no longer the full 
condensate fraction, but only its k = component. As an example, consider the 2D lattice 
potential (14) of figure 1, and compute the angle-averaged OBDM p(s) = (2n)^ 1 j d9p(s) 
[7]. The inhomogeneous mean-field contribution (15) then averages to \ ^ • V^Jo(Kjs), 
with Jq(.) the Bessel function. This function indeed goes to zero as s — > 00 and thus 
yields, by construction, the translation-invariant component as a result. But working with 
an configuration-averaged OBDM does not prove that the zero-momentum component is a 
good indicator for BEC, since the argument would be circular. Within mean-field theory, the 
condensate fraction is unity, regardless of the specific form of the condensate mode $(r). 
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The same argument applies whenever an ensemble-average over a random potential 
distribution is performed. In fact, the supposed disorder-depleted density calculated by 
Huang and Meng [14] and Giorgini, Pitaevskii, and Stringari [15] is nothing but the averaged 
condensate deformation n c V2 for a white-noise disorder potential. Within their approximation 
scheme, Huang and Meng were not able to calculate the condensate depletion induced by the 
external potential; this is achieved, for the first time to our knowledge, in Sec. 4 below. 

Summarizing the mean-field discussion, we reiterate that Bose-Einstein condensates 
form in many different spatial shapes and with different momentum distributions, determined 
by the competition of kinetic, interaction, and potential energy. In a trap or any spatially 
inhomogeneous potential, the zero-momentum eigenstate relevant to free space has no reason 
to determine the condensation into a single-particle orbital. In short, the population of 
momentum k = is not a good indicator for BEC in inhomogeneous systems. 

As a consequence, the population of momenta k ^ does not measure the depletion 
of inhomogeneous condensates, contrary to what has been suggested, unfortunately, in the 
groundbreaking work of Huang and Meng [ ], followed in this respect by numerous others 
[5, 15-21]. In order to assess the inhomogeneous condensate depletion, one requires a theory 
for the fluctuations. 

3. Momentum distribution of quantum fluctuations 

In this section, the momentum distribution Srik of fluctuations, as defined by (6), is determined 
via Bogoliubov theory as a function of the external potential. For a weak potential, a 
perturbative, but fully analytical, expression is obtained. 

3.1. Bogoliubov excitations of an inhomogeneous condensate 

Fluctuations around the inhomogeneous ground state <£>(r) are best described using the 
density-phase representation: (^(r) = cxp{iSip(r)}[^(r) 2 + Sfifr)] 1 ' 2 — $(r) develops 
as 

S$f(r) = —$(r)5h(r) + i$(r)5<p(r) + ... (16) 
2n c 

Here, the inverse condensate amplitude <t(r) — n c /$(r) is well defined because weak 
external potentials do not fragment the condensate. Likewise, highly excited states such as 
vortices are not considered, and so $(r) > holds everywhere. 

Fluctuations define Bogoliubov quasi-particles, or "bogolons". Mathematically, these 
are obtained by the canonical transformation [15] 

7 P = Sfip/ (2a p n 1 J 2 ) + ia p nl /2 Sip p . (17) 

Here, a p = (ep/ep) 1 ^ 2 for all p =/= is the traditional Bogoliubov transformation parameter, 
given by the ratio of free-particle dispersion e p — h 2 p 2 /2m to the Bogoliubov dispersion 
e P = K(4 + 2 9 n c)} 1/2 - For the zero mode, it is appropriate to define clq = 1, as discussed in 
Appendix A. Importantly, the fluctuations 5h p and 5tp p are the Fourier components of density 
and phase deviations away from the deformed mean-field ground state <£>(r) — and not from 

1/2 

the homogeneous background n c since then the deformation effect of the potential would 
be missed. By consequence of (17), the fluctuation (16) is expressed via bogolons as 

S^k = ^ (ukp-jp - Vkpj-pj , (18) 
p 
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where the inhomogeneous Bogoliubov transformation matrices 

Uk P = 2^fpf [°p I$fc -P + a P®k-p] > ( 19 ) 

Vkp = 2^fpf [ a P I$fc -P - a P*fc-p] > (2°) 

contain the Fourier coefficients of the condensate amplitude, = L~ d / 2 J d d re~ lk r Q(r), 
and its inverse, = [n c /<!>]*;. Some useful properties of this transformation to Bogoliubov 
quasiparticles are discussed in Appendix A. 

Inserting (18) and its Hermitian conjugate into Sn k = (6^ k S^fk) brings the fluctuation 
momentum distribution in the form 

Sn k = {"WKpI 2 + (u*kp u k P > + v* kp v kp ,) (7^V) - ( u kp v k P > (7p7- P ') + c - c -)} -(21) 

P:P' 

In principle, this equation is exact, and holds to arbitrary order in potential strength V. In order 
to compute the expectation values (j p j p >) and (7 P 7l p /), we need to specify the Hamiltonian 
of inhomogeneous Bogoliubov fluctuations. 

3.2. Inhomogeneous Bogoliubov Hamiltonian 

The quadratic Hamiltonian for the Bogoliubov excitations of an inhomogeneous Bose gas 
was derived in [9, 22] by a saddle-point expansion of the many -body Hamiltonian around the 
deformed ground-state solution $(r): 

H = J2 + J2 ffcVhfc/f w . (22) 

k k,k' 

The Bogoliubov-Nambu (BN) pseudo spinor = (jI , j-k) / V% allows for a rather compact 
notation. The only approximation in the derivation of the Hamiltonian is the neglect of third 
and fourth order terms in the fluctuations. In contrast, (22) is still exact in the external 
potential strength and has the structure H = + H^ v \ The price to be paid for spatial 
inhomogeneity is the appearance of the effective scattering vertex 

v ™ = ( Z Z ) ■ 

Its matrix elements, 

W k k' = t [akdk'Rkk' + a k ~ 1 a k ~, 1 Skk'] , (24) 



Ykk' — 7 [akdk'Rkk' — a k 1 a k , 1 Skk'] , 
are entirely determined by mean-field amplitudes and via 

Skk' = jd£~ k ■ k '( i - <w)X! $ fc-p $ p-fc' 

p 

25 



(25) 



(26) 



Rkk' = j^ d e ^2[k-k' + (k + k'- 2pf] * fc _p*p_ fc , - 2e° k 5 kk '. (27) 



p 



Here, we have dropped the superscripts and R^ v ^ used in [ ]. 
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3.3. Bogolon populations 

It is in principle possible to diagonalize the Hamiltonian (22) numerically, for each realization 
of the external potential, after having solved the nonlinear GP equation (7). However, for the 
purpose of analytical calculations in weak external potentials, a more economic strategy is to 
calculate the bogolon populations required in (21) perturbatively. We assume an equilibrium 
state at finite temperature T. The normal expectation value can be expressed via the single- 
quasiparticle Matsubara-Green (MG) function: 

ill'lp) = - G pp > (t) = --^ G pp , (iuj n ) (28) 

where w„ = 2irn/(3 are the bosonic Matsubara frequencies related to the inverse temperature 
f3 = l/k-oT. Similarly, the anomalous expectation value 

is expressed in terms of the anomalous MG function F(z). Together, the normal and 

anomalous MG function enter the Nambu-MG matrix Q = ( _ G t ) which expands as 

G = Go + GoVGo + GoVGoVGo + ■ ■ ■ The free propagator is diagonal in Nambu space 
and momentum representation, Go P {iu n ) — diag(Go p («w„), Gq p (— iuj n )) with Gq p {z) = 
[z — e p ] _1 . Matsubara sums like (28) and (29) are carried out using textbook recipes such 
as (1 1.58) in [23]; each simple pole with energy e p contributes one Bose-Einstein occupation 
number 

v := K/3e P ) = [exp(/3e p ) - l]" 1 . (30) 

The expectation values (28), (29) are then straightforward to calculate. For brevity, we 
present here only the diagonal result for p' = p, i.e. the bogolon population v p = (7 P 7 P ), up 
to order V 2 : 



de ^— ' I e — e' \de 

p' 

1 (dv 1 



VP 



Ypp'Yp'p ( • (31) 



e + e' \ de e' + e 

Here the short-hand notations e = e p , e' = e p >, and v' — v(/3e') are used. At temperature 
T = when all occupation numbers and their derivatives vanish, there is only a single finite 
contribution to the normal bogolon population due to the external potential, to order V 2 : 

V P = u T e r\2 PP P'P- @2) 

p , 

The bogolon quasi-particles are populated by the random potential even at zero temperature 
because the full Hamiltonian (22) is not diagonal in the basis that diagonalizes . 
Similarly, the anomalous population ir p = — ( A fpJ- p ) reads, to order V 2 , 

_ 1 + 2. _ 1 ^ e - e 1 + 2(a/ - eV) 
n P 2e pp e ^ (e-e>)(e + e>) Y pp' W p'p- (ii > 

At T = 0, the result simplifies slightly: 

n p = 2e^PP — g $Z e _)_ g7 Ypp'Wpi p . (34) 
p' 
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With this, everything is in place to calculate the full momentum distribution (21), or 
equivalently the OBDM (5). In the following, we pursue a fully analytical calculation by a 
perturbative expansion up to order V 2 in the bare external potential. In the upcoming section 
3.4 the Bogoliubov transformation matrices (19) and (20) are determined together with the 
scattering matrix elements (24) and (25). In Sec. 3.5 the results are collected into a compact 
expression for the momentum distribution. 



3.4. Weak-potential expansion 

The perturbation expansion (8)-(9) of the condensate amplitude Q k in powers of V implies a 
similar expansion for the Bogoliubov transformation matrices (19) and (20): 

«fcp = «S+«S +u S+-" (35) 

«*p = < } +«£+«$ + ••• C 36 ) 
To zeroth order in the external potential, the transformation matrices are diagonal in 
momentum, as required by translation invariance, and the traditional Bogoliubov amplitudes 
are recovered: 

1 



v 



(0) 



-(a p 1 + a p )S kp = u p S kpi (37) 



v kl = \( a p 1 _ a p )8 kp = v p 5 kp . (38) 
To first order, the matrix elements are proportional to the potential matrix element (10): 

u kp = -v P Vk-p, 4p = -UpV k - p . (39) 
For the second-order matrices, only the diagonal matrix elements will be required: 

2 

where V<z of Eq. (12) is second order in V. 

As spelled out in [ ], also the BN scattering potential (23) admits the expansion 
V = V^ 1 ) + V^ 2 ) + . . .. Also here, the first-order scattering amplitudes 

Wff = gn c w^V k . p , = gn c y^V k . p (41) 

are proportional to (10), with 



U kk = o V *> Vfcfc = = V 2 , (40) 



= ? [o-ka p {k 2 +p 2 -k-p)-a k x a p [ k • p] 



= £ 2 [(uk - v k )(u p - v p )(k 2 +p 2 ) - 2(u k u p + v k v p )k ■ p] , (42) 
Vki = ¥ [ a kap(k 2 +p 2 -k-p)+ a-'a^k ■ p] 

= C 2 [i u k - v k ){u p - v p )(k 2 +p 2 ) + 2(u k v p + v k u p )k ■ p] . (43) 
Second-order scattering amplitudes are later only needed for k = p: 

W^^Yf^gn^&H-v? (44) 
p 

with 

tfi =2a 2 e[k 2 + (k-p) 2 ]. (45) 
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3.5. Momentum distribution 

Collecting all results up to order V 2 , the single-particle fluctuation momentum distribution 
(21) reads 

5n k = (5& h 6* k ) = 5nf + 8nf , (46) 

where the superscript indicates the order in the external potential strength V. To zeroth order, 
i.e. for the homogeneous system, we recover the well known zero-temperature momentum 
distribution [8] 

as a consequence of the two-body contact interaction. The first-order term vanishes by 
momentum conservation, Sn^ — 0. The external potential induces the second-order shift 

*n<?=£Mg|V fc _ p | a , (48) 
p 

whose kernel is defined in terms of (42), (43), and (45): 



= (vl - 2u k v k + Up) - 2(u k u p + v k v p )- — — 



(k + ( P 



+ f..a , v 2 ^n c y k p) gn c u k v k f (2) 2gn c (1) (1) 1 

This expression, together with the preceding general form (21), constitutes the main result of 
this calculation. From here on, we explore its consequences by studying two generic examples 
of external potentials: a weak lattice potential on the one hand, and a random potential on the 
other. 



5.5.7. Lattice potential - A pure lattice potential like (14) has only the Fourier components 

v k = \ E?=i v j ( S k,Kj + &h-K s ), such that 



1 d 



\ Vk \ 2 = i V J ( Sk ^ + S K-K t ) ■ (50) 

3=1 

Thus, the momentum distribution shift (48) is given by 



5 



k 



(2) _ 1 (Vj/gn c ) 2 / - ( 2) ~ (2) \ 

A \0 + (TC.W\2 l-^fefe+JCf + M hk-K 4 ■ W 



4 [2 + (if^) 2 



Figure 2 shows the momentum distribution (46) in the square lattice potential of Fig. 1 . 
The top panel shows a cut through the distribution at k y = 0, together with the separate 
contributions of the homogeneous and potential-induced distribution. Compared to the mean- 
field figure 1, the quantum fluctuations broaden the momentum distribution substantially. 

In formula (5 1 ), the product Kj£ compares the characteristic length scale of the potential, 
KJ , with the condensate healing length £. For not-too-low densities and typical interaction 
strengths achievable with ultracold atoms, one easily reaches Kj£ <C 1, known as the 
Thomas-Fermi (TF) regime. The healing length is the characteristic scale also for the entire 
kernel (49). In the deep TF regime Kj£ — > 0, and for finite momenta k > Kj, this 




Figure 2. Zero-temperature momentum distribution <5nj. of quantum fluctuations, Eq. (46), 
in the square lattice potential of Fig. 1. The upper panel shows a cut along k y = 0. The 

main contribution stems from the isotropic background Sn^ , Eq. (47). The potential-induced 

distribution Sn k , Eq. (51), reflects the lattice structure. 



~ (2) 

complicated kernel can be approximated by the diagonal term . The potential-induced 
change of momentum distribution then takes the simple isotropic form 

Sn^ - V - (fc ^" 1 (52) 

Here, v 2 — ^2~2j(Vj I ff n c) 2 measures the potential variance in units of the mean-field 
interaction energy. In the TF regime (52), the external potential is found to shift population 
from low momenta fc£ < 1 to high momenta fc£ > 1, and this independently of the detailed 
form of the potential. 



3.5.2. Random potential - A random potential can be seen as a superposition of many 
lattices with a random distribution of Fourier components V^, specified by the ensemble 
averages Vk, VkV p , etc. Here, we assume without loss of generality that the potential is 
centred, V(r) — or Vb = 0. All we need at order V 2 then is the pair correlator 

VqVlq, = L-% q ,V 2 a d C d (qcr). (53) 

The dimensionless function Cd(qo-) characterizes the potential correlation on the length scale 
er; the normalization is chosen such that (a/L) d 2~2 q Cd{qcr) — 1 in the thermodynamic limit. 
Using (53), the ensemble-averaged change of the single-particle momentum distribution (48) 
takes the form 




Figure 3. Quantum-fluctuation momentum distribution shift (54) induced by a 2D random 
potential of reduced variance v 2 = V 2 / (gn c ) 2 and Gaussian correlation, (55), for different 
values of the correlation length relative to the condensate healing length, £ = <x/£. For 
comparison, also the momentum distribution Sn^ = v 2 , Eq. (47), of the homogeneous 
fluctuations is shown. The insets show the normal and anomalous populations of Bogoliubov 
excitations, (56) and (57), respectively, also as function of k£. In the TF limit £ — > oo, 
all distributions converge towards the universal expressions (52), (58), and (59), respectively, 
shown in dashed black. 



where v 2 = (V/gn c ) 2 is the potential variance in units of mean-field interaction energy. 

Figure 3 shows the ensemble-averaged, isotropic, momentum distribution (54) that is 
induced by a random potential with Gaussian correlation 

C d [qa) = (2Ti) d ' 2 exp{-g V/2} (55) 

in d = 2 dimensions as function of |fc|£. Different curves correspond to different values of 
the correlation parameter £ = cr/£, namely the correlation length relative to the condensate 
healing length. Just as for the lattice, also here things simplify considerably in the TF regime 
<J 3> £ where the disorder correlation length is much longer than the condensate healing 
length. Then, the potential correlator tends to a ^-distribution, and (54) reduces to (52), plotted 
in dashed black. Clearly, the momentum distribution of quantum fluctuations is given by the 
universal form (52) in any external potential that is sufficiently smooth to yield a Thomas- 
Fermi condensate profile. 

The insets of Fig. 3 show the ensemble-averaged, normal and anomalous bogolon 
populations i>h = {iklk) and ir k = — {iki-k) at zero temperature. These populations, 
(32) and (34), are given by 



*fc=^ 2 T7> q , 2 m , (56) 




(yifc-g) 2 C d (q<7) 
q j [2 + (<?0 5 



r-— \ > yik-q w kk-g \ C d {qa) 

' } S 2gn c e k e k (e k + e fc _ q ) f [2 + (q^) 2 } 2 ' K ' 
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in terms of the envelopes (42), (43), and (45). In the TF regime, they tend toward the universal 
limiting expressions (dashed black in the insets of Fig. 3) 

_, 2 [l + (fcQ 2 ] 2 

^ TF ~T[2+(fcOT' ( } 



v 2 2 + 4(fcQ 2 + (fcg) 4 

TTfcTF = —, ro i n ,-\9i9 ' ^> 



4 [2 + 

and this independently of the potential details. 



4. Quantum depletion of the condensate 

The total particle density p(0) = n = n c + Sn is the sum of condensate density n c and 
non-condensed density Sn. The condensate fraction is n c /n = 1 — Sn/n. Within Bogoliubov 
theory, the existence of a finite non-condensed fraction Sn/n at temperature T = is called 
"quantum depletion", because it arises from quantum fluctuations around the mean-field 
approximation to the true condensate. From a many-body point of view, the non-condensed 
fraction is of course not more quantum than the condensed one, or perhaps even rather less. 
Here, we follow the established nomenclature and continue to speak of quantum depletion, at 
zero temperature, as opposed to the thermal depletion at finite temperature. 

From definitions (5) and (6) it follows that the depleted density Sn — n — n c is the 
integral of the fluctuation momentum distribution, 

Sn = L- d ^5n k . (60) 

k 



4.1. Homogeneous system 

Let us first recall the homogeneous case V = [ ]. Since condensation occurs in the k = 
mode, the depleted density simply contains all particles with finite momenta, and the zero- 
temperature momentum distribution (47) implies 

6n<® =L- d Y,vl (61) 

In d = 3, Eq. (61) evaluates to the depleted density Sn^ = [eV^^ 2 ^ 3 ]" 1 in the 
thermodynamic limit. Equivalently, the relative depletion reads Sn^ ' jn — 8(na 3 ) 1 / 2 /37r 1 / 2 
because g = 4irh 2 a s /m in terms of the s-wave scattering length a s and £ 2 = ti 2 / (2mgn) (to 
leading order, we can identify n c w n in all perturbative results). The Bogoliubov ansatz is 
justified whenever the fractional depletion is small, 5n(°> <C n or equivalently n£ d ^> 1. This 
is the case whenever the so-called gas parameter na? s is small, i.e., for low enough density or 
weak scattering. 

In d = 2, one finds Sn^ = [Stt^ 2 ]^ 1 . The quantum depletion Sn^ jn is roughly 
independent of density, and requires weak scattering. 

In d = 1, the infrared fc _1 -divergence of Vk under the integral prevents the existence of 
a homogeneous ID condensate. But quasi-BEC occurs in finite geometries. Cutting off the 
ID-integral at some value a = £fciR <§C 1, with km of the order of the inverse system size, one 
finds <5n(°) = (2 In 2 — 2 — \na)/(2y/2ir{;), up to order a. Bogoliubov theory then is valid 
whenever ^ 1, i.e., requires high enough density in order for the mean-field picture to 
apply in the first place (in the very low-density limit one crosses over to the Tonks -Girardeau 
regime). 
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So in all relevant dimensionalities, there is a window of validity for Bogoliubov theory, 
and the depleted density can be written 5n^ ' = c c [£~' i , with a d-dependent numerical constant 
Cd of order unity or smaller. 

4.2. Potential depletion 

In an inhomogeneous system, the quantum depletion cannot be calculated by counting all 
particles with finite momentum, as argued in Sec. (2) above. Instead, the depleted density 
(60) is the integral of the fluctuation momentum distribution (21), which splits into two 
contributions: the quantum depletion of the homogeneous system plus the potential-induced 
depletion properly speaking. We like to emphasize once more that the reason for the potential- 
induced depletion is not simply the appearance of k ^ components in the momentum 
distribution. In traps, lattices, or random potentials, already the condensate contains finite 
momentum components, Eq. (11). And clearly, the condensate fraction of an inhomogeneous, 
but non-interacting BEC (such as noninteracting cold atoms in a harmonic trap) is unity, 
regardless of its spatial configuration. Rather, an external potential can change the condensate 
fraction because it modifies the local particle density. This change in particle density changes 
the local interaction energy, which in turn changes the depletion. Since the interacting system 
has a nonlinear response, even in a purely sinusoidal lattice potential the high-density regions 
will deplete more condensate than the low-density regions can gain back. At the end, the 
presence of the potential causes a net additional depletion of the condensate, an effect that we 
propose to call "potential depletion". 

To our knowledge, the potential depletion, beyond the mean-field deformation of the 
condensate, has never been calculated analytically. In approaches very similar to ours, Singh 
and Rokhsar [24] arrived at numerical results for the potential depletion; Lee and Gunn [25] 
estimated a different depletion. Within our inhomogeneous Bogoliubov theory, computing 
the potential depletion is straightforward: Using the perturbative result (48) in (60), we find 
for the potential-depleted density 



The sum over k can be carried out without touching the potential, whence the second equality, 
which defines the isotropic depletion kernel for the bare potential, 



Prefactors are chosen such that in the thermodynamic limit Md(q£) is a dimensionless 
function of q£ only. 

4.2.1. Lattice potential - For the lattice potential (50), the potential depletion (62) reads 



For the 2D lattice potential of Figs. 1 and 2, one finds Sn ( ^ = M 2 {K£)V 2 / '(£#n c ) 2 m 
0.141 Sn^ '; for these parameters, the potential depletion amounts to only 14% of the 
homogeneous depletion. 

Let us furthermore check (64) against the QMC results of Astrakharchik and Krutitsky 
[7], who investigated two different interaction strengths in a square lattice, such that K£ — 





(63) 




(64) 
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Figure 4. Disorder-induced quantum depletion (67), relative to the clean value and in units of 
disorder strength v' 2 , as function of the correlation ratio f = for the Gaussian correlation 
(55). The curve for d = 1 depends weakly on the infrared cutoff a = £&ir that regularizes 
already the clean case; a = 0.05 in this plot. Data points are obtained by an exact numerical 
diagonalization of the Bogoliubov Hamiltonian (22), followed by an ensemble average over 
disorder. Error bars denote the estimated error for the average; the number M of realizations 
is chosen for each point such that ML/ a = 2000. Dashed: universal Thomas-Fermi limit 
Atf =3/8 for very smooth disorder (f — > oo) in d = 3. The curves for lower dimensions 
tend to Atp = and Atf = —1/8, in d = 2 and d = 1, respectively. Dotted: limiting 
behaviour for 5-correlated disorder (C < 1): A(£) = t3 d ( d C d {0) with ft « 0.245 (for 
a = 0.05), #2 ^ 0.135, /3 3 » 0.160. 



6.275, and = 1.984, respectively. For these values, (64) predicts a potential depletion 
of 5n^ = 1.05 6nl°\ and 0.35^°), respectively. If one takes the QMC values for 
the homogeneous depletion as the reference, then the final condensate fraction should be 
Nq/N = 1 — Sn/n = 0.98 in the one case and 0.73 in the other, which is in good agreement 
with the data [ ], No /N ss 0.99 and 0.7, respectively. Incidentally, the Bogoliubov prediction 
for the clean depletion Sn^ / n does not agree so well with the data, which may be due to finite 
size effects or the slightly different interaction potential used in the QMC approach (hard-core 
bosons instead of s-wave scattering). 



4.2.2. Random potential - Using the momentum distribution (54) in (60), the quantum 
depletion (62) by a random potential with correlation (53) is found to be 

^ (2) = JCwl* E MMMqa). (65) 

Scaled by the homogeneous depletion in the thermodynamic limit, this can be written as 

SnW = v 2 Sn^ A(C) (66) 
where v — V/ (gn c ) and 

C d f d d u 

A ^ ^ ~ / 7^- d Md(u)C d (u(). (67) 
c d J (2tt) !1 

This relative potential depletion is found to be a function of the correlation ratio Q = <r/£. 
Only in d = 1, it depends also very weakly on the cutoff a that regularizes already the 
clean depletion. The convergent integrals (66) require no additional ad-hoc cutoffs, neither 



Bose-Einstein condensates in external potentials 



15 



infrared (since the excitations are orthogonal to the vacuum) nor ultraviolet (since potential 
correlations are included). Fig. 4 shows A(£) for a Gaussian-correlated random potential 
(55) in dimensions d = 1,2,3. Data points result from the numerical diagonalization of 
the Bogoliubov Hamiltonian (22) in a system of linear size L such that a = £/L = 0.05, 
followed by an ensemble average over disorder. The exact shape of the curve depends on the 
correlation function, the general features, however, are rather robust. 

The asymptotic behaviour of A(£) for very small or very large correlation lengths is 
simple. In the (5-correlated limit £ — > of a white-noise potential, the generic scaling of 
(67) is A(C) = f3 d ( d C d (0) with B d = c^ 1 J -0fM d (u). The numerical coefficients are 
Pi R3 0.245 (weakly dependent on the cutoff a), /?2 ~ 0.135, /3s 0.160. In this white-noise 
regime, the depletion depends on a d C d (0), and thus requires the existence of a microscopic 
correlation scale. In the opposite limit £ — > oo of the Thomas-Fermi regime, the result 
converges to a truly universal limit 

A TF = c-M d (0) = ^/ ( g,^ F (68) 

and evaluates, using (52), to Atf = 3/8 in d — 3 and Atf = in d — 2. In 2 and 3 
dimensions, the depletion is non-negative, as one would expect for a random potential that 
should broaden the momentum distribution overall. For d — 1 the TF-limit evaluates to the 
negative value Atf = —1/8 (in the limit of infinite system size). This would seem to imply 
that the random potential re-populates the condensate. But also in d = 1 the depletion is 
positive for most values of (, as shown in Fig. 4. The curve only crosses over to negative 
values for such a large value £ = er/£ (depending on the cutoff a), that the correlation length 
a has to be comparable to the system size, which is not the regime of present interest. 

Interestingly, the TF limit for the potential depletion can also be derived by the local- 
density approximation (LDA) i%tf =n c — V(r)/g combined with the scaling 6n = c d £~ d = 
c ' d n ^ 2 °f tne homogeneous depletion (Sec. 4.1): 

Sn^F = c' d n d J 2 [l-V(r)/gn c } d / 2 = Sn^ 

and thus Atf = d(d — 2)/8, in agreement with the result of (68). This argument shows that 
in d = 2 dimensions, the TF potential depletion is zero even non-perturbatively since V = 
without loss of generality. This LDA reasoning works for the correction of the depletion, 
but not for the excitation dispersion relation, where genuine scattering effects determine 
corrections to the speed of sound and density of states [9, 26]. 

Summarizing the results of this section, we conclude that the combined depletion due to 
interaction and external potential reads 

Sn = Sn (Q) [l + v 2 A], (69) 

with |A| <C 1. Clearly, the potential depletion alone, SrS 2 ^ jn — (5n^ / n)v 2 A, is at least 
a factor Sn^/n -C 1 smaller than the mean-field condensate deformation (11), which is 
of order v 2 . In hindsight, this result is rather plausible: the primary effect of the external 
potential is merely to deform the condensate. The potential depletion is a secondary effect, 
caused by enhanced interaction in the regions of higher density. We conclude that, as long 
as the original assumption of a non-zero condensate amplitude holds, our inhomogeneous 
Bogoliubov theory applies to Bose condensates in rather inhomogeneous potentials. 



1 d{d-2) 2 _ 3 , 
1 + — '-v 2 + 0(v 3 ) 
8 



Bose-Einstein condensates in external potentials 



16 



5. Summary 

We have investigated the effect of external potentials on Bose-condensed gases using 
inhomogeneous Bogoliubov theory. The principal effect of an external potential is to deform 
the mean-field condensate. Secondly, the potential affects the momentum distribution of 
quantum fluctuations, for which we have obtained a general expression. Finally, we have 
calculated the quantum depletion induced by the external potential, or potential depletion 
for short. In detail, we have studied lattices and spatially correlated random potentials. 
The potential depletion turns out to be proportional to the homogeneous depletion, a 
fact that underscores the applicability of inhomogeneous Bogoliubov theory in weak to 
moderately strong potentials. Our analytical predictions are in agreement with a numerical 
diagonalization of the Bogoliubov Hamiltonian as well as with recent quantum Monte Carlo 
simulations [ ]. The inhomogeneous Bogoliubov theory shown at work here is therefore 
proven capable of describing the excitations of weakly interacting condensates in external 
potentials, and from there ought to provide many of other static and dynamic properties. 
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Appendix A. Transformation to Bogoliubov quasi-particles 

The matrices u kp and v kp as defined by (19) and (20) are the Fourier components with wave 
vector k of the modes u p (r) and v p (r) defined in [9]. As explained there, the momentum 
index p can be used to label the modes even in the inhomogeneous setting. The transformation 
(18) preserves the canonical commutation relation, and thus guarantees [7 P ,7 P /] = S PP ' as 
well as [j p , j p i] — 0, via the completeness relations 

^2 ( U kp u *k'p - v kpV* k , p ) = S k k', (A.l) 

p 

( U kp V k'p - V kpU*k'p) = °- ( A - 2 ) 

p 

The non-symmetric matrices Uk P ^ u pk and v kp ^ v pk also satisfy the biorthogonality 

^2(u kp u* kp , - v kp vl pl ) = 5 pp ,, (A3) 

k 

J2(u kp v* kp , - v kp u* kpl ) = 0. (A.4) 

k 

The zero mode deserves special attention because a p = (e^/cp) 1 / 2 diverges as jT 1 when 
p^O. In this range, elementary excitations are essentially phase fluctuations. Setting ao = 1, 
one finds that the Bogoliubov excitation 

70 = 5n /{2n 1 J 2 ) + in^Sipo, (A.5) 

together with its Hermitian conjugate, describes the number fluctuation 



SfiQ = ■y/n^(% +7o) = L i l d d r$(r) <5v&(r) + 5&{r) 



(A.6) 
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This operator, called P in [ ], generates an exact zero-energy (Goldstone) mode of the U(l)- 
symmetry breaking Bose condensed state. The corresponding mode functions are 

u ka = [$ fe + $ k ]/{2N^ 2 ), v k o = [$ fc - $ fc ]/(27V c 1 / 2 ). (A.7) 

With these definitions, the completeness relations (A.1)-(A.2) and biorthogonality relations 
(A.3)-(A.4) include and extend to the zero modes. In the present article, we investigate the 
spatial structure of quantum fluctuations, and the contribution from p = has vanishing 
weight anyway in the thermodynamic limit where sums over momenta turn into integrals — 
except for ID, but there, we introduce an IR cutoff to regularize the divergence. This masks 
the phase diffusion physics at long distances and times [8, 27, 28], which has not been the 
subject of the present investigation. 

References 

[1] M. Lewenstein, et aL.Adv. Phys. 56, 243 (2007). 

[2] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 

[3] L. Sanchez-Palencia and M. Lewenstein, Nat. Phys. 6, 87 (2010). 

[4] O. Penrose and L. Onsager, Phys, Rev. 104, 576 (1956). 

[5] G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, Phys. Rev. A 66, 023603 (2002). 

[6] S. Pilati, S. Giorgini, andN. Prokof'ev, Phys. Rev. Lett. 102, 150402 (2009). 

[7] G. E. Astrakharchik and K. V. Krutitsky, Phys. Rev. A 84, 031604 (201 1). 

[8] L. Pitaevskii and S. Stringari, Bose-Einstein condensation, Oxford University Press, New York (2003). 

[9] C. Gaul and C. A. Muller, Phys. Rev. A 83, 063629 (201 1). 

[10] F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 (1996). 

[11] L. Sanchez-Palencia, Phys. Rev. A 74, 053625 (2006). 

[12] P. Lugan and L. Sanchez-Palencia, Phys. Rev. A 84, 013612 (201 1). 

[13] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. Bloch, Nature 415, 39 (2002). 

[14] K. Huang and H.-F. Meng, Phys. Rev. Lett. 69, 644 (1992). 

[15] S. Giorgini, L. Pitaevskii, and S. Stringari, Phys. Rev. B 49, 12938 (1994). 

[16] M. Kobayashi and M. Tsubota, Phys. Rev. B 66, 174516 (2002). 

[17] G. M. Falco, A. Pelster, and R. Graham, Phys. Rev. A 75, 063619 (2007). 

[18] V. I. Yukalov and R. Graham, Phys. Rev. A 75, 023619 (2007). 

[19] Y. Hu, Z. Liang, and B. Hu, Phys. Rev. A 80, 043629 (2009). 

[20] S. Pilati, S. Giorgini, M. Modugno, and N. Prokof'ev, New Journal of Physics 12, 073003 (2010). 

[21] C. Krumnow and A. Pelster, Phys. Rev. A 84, 021608 (201 1). 

[22] C. Gaul and C. A. Muller, Europhys. Lett. 83, 10006 (2008). 

[23] H. Bruus and K. Flensberg, Many-body quantum theory in condensed matter physics, Oxford Univ. Press 
(2004). 

[24] K. G. Singh and D. S. Rokhsar, Phys. Rev. B 49, 9013 (1994). 

[25] D. K. K. Lee and J. M. F. Gunn, /. Phys.: Condens. Matter 2, 7753 (1990). 

[26] C. Gaul, N. Renner, and C. A. Muller, Phys. Rev. A 80, 053620 (2009). 

[27] M. Lewenstein and L. You, Phys. Rev. Lett. 77, 3489 (1996). 

[28] L. Fontanesi, M. Wouters, and V. Savona, Phys. Rev. A 81, 053603 (2010). 



